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Due to entire pneumatic connection between free piston and free displacer, the motion parameters of 
them including amplitude and phase shift can actually impact the cooling capacity and overall perfor¬ 
mance of cryocooler obviously. In this study, the procedure of design and manufacture pneumatic free 
piston and free displacer (FPFD) Stirling cryocooler had firstly been described in details. Then in order 
to accomplish study, the experimental bench has been set up based on 80 K@1 W Stirling cryocooler. 
The effect of the thermodynamic and pneumatic parameters including charging pressure, natural fre¬ 
quency of displacer, damping coefficient of displacer, working frequency on the pressure, displacement 
and displacer phase shift has been investigated, respectively by means of experimental and theoretical 
method. In particular, the variation of damping is realized by adjusting the width of clearance cut on 
the additional damping component, which is screwed on the displacer rod. Similarly, natural frequency 
of displacer is changed by the extra mass connected on the displacer. Due to the results of experimental 
study, the optimum working conditions of this Stirling cryocooler for 80 K cold tip temperature are as fol¬ 
lows: charge pressure 15 bar, natural frequency of displacer 46 Hz, width of clearance 300 pm and work¬ 
ing frequency 43 Hz. In agreement with the optimum working conditions, neighborhood interval of 90° is 
the ideal working domain for displacement phase shift. Meanwhile, the displacer phase shift should 
approach to 0°as near as possible and pressure phase shift should also be as small as possible, which have 
linear relation with non-dimensional damping characteristic of compressor. In view of theoretical study, 
the expressions of three phase shifts deduced from thermodynamic equation of piston and displacer 
respectively are expressed as the functions of working parameters, which are verified by the experimen¬ 
tal data and consequently can be used as the powerful guidance to optimum seeking. 

© 2008 Elsevier Ltd. All rights reserved. 


1. Introduction 

The growing demand for long-wavelength infrared imaging 
instruments for space observational applications, high tempera¬ 
ture superconducting filters and military detector has led to the 
ongoing development of several second generation Stirling cycle 
cryocoolers that are scaled up in size from the original Oxford de¬ 
sign. Recently, fast cool down, small size, lightweight, low power 
consumption and high reliability split-type Stirling cryocoolers dri¬ 
ven by linear compressor have been used to meet the long operat¬ 
ing lift requirements. The free piston free displacer (FPFD) Stirling 
cryocooler developed initially by the Philips laboratories consists 
of two-compressor piston driven by linear motors, which make 
pressure waves and a pneumatically driven displacer piston with 
a regenerator. In general, the efficiency of the Stirling cryocooler 
is mainly affected by the efficiency of the linear motor, the reso¬ 
nant frequency of the compressor and expander, the displacement 
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of the piston and the displacer, and the phase shift between the 
piston and the displacer. 

As shown in Fig. 1, in Stirling cryocooler, the displacer in the 
expansion cylinder is free and oscillates simply as a result of the 
gas force acting upon it. This gas force is the direct result of motion 
of compressor piston driven by the electric force. Due to the com¬ 
plex characteristic of gas spring in the warm and compression 
space, the phase shift between the displacement of piston and dis¬ 
charge pressure can reveal the performance of linear compressor to 
some extent. Moreover, the gas force provided by the actuating 
piston oscillating in the compression cylinder with a phase shift 
relative to the displacement of displacer indicates the volume 
change in the expansion space, and its pressure-volume diagram 
illustrates the work done by the gas on the displacer. This positive 
work result in the refrigeration effect and the cold finger becomes 
cold. Finally, the displacement phase shift between piston and dis¬ 
placer reflects the relative position of two moving components lit¬ 
erally, but it actually has great effect on the cooling capacity and 
electric power consumption of cryocooler. 

These three phase shift characteristics indicate performance of 
linear compressor and cryocooler from different viewpoints. Park 
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Nomenclature 



A 

area (m 2 ) 

(P 

phase shift between expansion pressure and displacer 

B 

magnetic density (T) 


(°) 

C 

damping coefficient (N s/m) 

y 

phase shift between the force and electric current (°) 

F 

amplitude of force (N) 

Subscripts 

/ 

harmonic electric force (N) 

b 

bounce space 

i 

current (A) 

c 

compressor space 

k 

stiffness (N/s) 

d 

displacer 

M, m 

mass (kg) 

e 

expansion space 

P 

pressure (Pa) 

ee 

electric 

Q 

refrigeration capacity (W) 

g 

gas spring 

t 

time (s) 

h 

hot space 

COP 

coefficient of performance 

m 

mechanical stiffness or damping 

X 

displacement (m) 

mean 

mean value 

X 

amplitude of displacement (m) 

ideal 

ideal value 

V 

buffer volume (m 3 ) 

n 

natural frequency 

Greek letters 

0 

mean value or operating frequency 

a 

poly tropic coefficient, 1 (isothermal) ^ a ^ y (adiabatic) 

P 

compressor piston 

CO 

angular frequency (rad/s) 

r 

rod 

-A 

displacement phase angle (°) 




[1] had measured the cool-down characteristics of the cold end 
with laser displacement sensor in the expander of the Stirling cryo- 
cooler. The amplitude of displacer’s pressure and displacement had 
been presented with the variation of charging pressure (15-25 bar) 
and operating frequency (40-55 Hz), when the displacement 
phases shift is 45°, optimum operating frequency was finally 
found, which was decided by many factors including resonant fre¬ 
quency of expander and friction and flow resistance. Kornhauser 

[2] had developed an improved linear dynamic model (anelastic 
model) of gas springs in the Stirling engine, which dealt the tradi¬ 
tional isothermal model in parallel with adiabatic model and could 
predict the phase shift relative to discharge pressure of compressor 
and piston displacement, stroke amplitude and damping character¬ 
istic theoretically. Gaunekar [3] had set up a mathematical model 
aiming to analyze the performance characteristic of cryocooler sys- 
temically in theory, although many assumptions had been 
adopted, some results can still be used as the guide of design. 

In this study, pressure, displacement and displacer phase shift 
characteristics have been investigated respectively by means of 
experiment under different conditions, such as working frequency, 
charged pressure, natural frequency and motion damping coeffi¬ 
cient of displacer. In addition, the corresponding theories of them 


have been deduced from dynamic equations of piston and displacer, 
which reveal the relationship between the working parameters and 
three phase shifts respectively. The experimental results of the per¬ 
formance and optimum operating conditions of the Stirling cryoco¬ 
oler according to the variations of the above phase shifts are 
presented as well. This phase shift is the powerful tool to the perfor¬ 
mance evaluation and optimum design of Stirling cryocooler. 

D i ~ p 1 Phase shift between piston displacement and discharge 
pressure (pressure phase shift \jj w ) 

D i ~ D 2 Displacement phase shift between piston and displacer 
(displacement phase shift i// pd ) 

Pi ~ D 2 Phase shift between discharge pressure and displacer 
displacement (displacer phase shift \j/ dd ) 

2. Design and manufacturing of the Stirling cryocooler 

2.1. Design of Stirling cryocooler 

In terms of the physical process, the pneumatic Stirling cryoco¬ 
oler can be regarded as the two-freedom damping vibration system 
as a whole. The working gas in compression, expansion chamber 
and regenerator accomplishes the thermodynamic process that is 



Fig. 1 . Schematic diagram of FPFD Stirling cryocooler. 
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actually coupling tightly with the dynamic working conditions. The 
equation describing the motion of piston and displacer [4] is as 
follows 


MfX(\ — (A e ^r)Ph A ePe "T ArPb ^md^d Qnd^d (1) 

MpXp = (Pb — p c )Ap — k m pXp — C mp x p + Bil (2) 

In above equations, the gas pressure, damping factor is essentially 
nonlinear term which needs very complex iteration process before 
getting the convergent solution. Getting precise mathematic solu¬ 
tion of whole Stirling cryocooler control equations is a very heavy 
burden to an engineer. In practical design, the isothermal model 
[5] always adopted to simplify the thermodynamic model and sub¬ 
sequently explicit expression of thermodynamic parameters can be 
obtained. Moreover, the motion of piston and displacer coupled 
tightly with each other by the pneumatic pressure term is always 
uncoupled by means of sine wave assumption of these two dis¬ 
placements. Meanwhile the first three terms in the right side of 
Eq. (1) and the first term in the right side of Eq. (2) reflecting the 
pneumatic force of expansion and compression space can be equiv¬ 
alently simplified to gas spring and equivalent hysteresis loss. The 
gas spring constant is given by [6] 



dP 


mean 


A 2 x 


A(x m - x) + V 



Then resonant frequency of piston and displacer is as follows: 

las'll < 4 > 

For the compressor, the gas stiffness is bigger than that of 
mechanical spring. So Eq. (3) should be involved in the natural fre¬ 
quency calculation according to Eq. (4). In contrast, the gas spring 
stiffness of displacer is weaker than that of spring, so it can be 
ignored. 

The ideal refrigeration (PV power) of cryocooler can be written 
as the following equation [1]: 

Qideai = nfP A AiX A sin cp (5) 


As shown in Eq. (5) the ideal refrigeration is greatly affected by the 
operating frequency, the amplitude of pressure, the amplitude of 


the displacer and the phase shift between the pressure and 
displacement. 

2.2. Manufacture of Stirling cryocooler 

The FPFD Stirling cryocooler for case study consisted of three 
major parts (Fig. 2): the linear compressor module, the expander 
module and electric power supply and control module. The linear 
compressor consisted of aluminum allay base, inner and outer iron 
yokes, permanent magnets (NdFe35), coils (copper wire), cylin¬ 
ders, pistons and flexure springs. The compressor was based on 
the well-proven ‘Oxford’ principles of spiral flexure springs and 
non-contacting clearance seals. The machine was a compact mov¬ 
ing piston design, with the piston and cylinder located within the 
core of the magnetic circuit of a moving coil motor. The springs 
were the only component subjected to significant fatigue loading, 
and there were routinely batch and sample tested at their full 
stroke in excess of 4 x 10 9 cycles. In addition, strain-stress test 
conducted indicated the linear stiffness relationship. The structure 
of the compressor consists of two opposite parts that were inher¬ 
ently well balanced with the two halves mounted in line and oper¬ 
ated in anti-phase. As a result, very low levels of self-induced 
vibration had been obtained. In terms of material selection, non- 
metallic materials used in the assembly such as plastics, adhesives 
and gasket were selected from an existing knowledge base of 
materials with a low out-gassing rate. As for metal material, the 
strength, stiffness, thermal expansion and processing characters 
had been considered and finite element analysis of some important 
components had even been done. Meanwhile, geometric toler¬ 
ances and surface requirement should commensurate with the 
functional demands [7]. 

The expander module consisted of a displacer, a regenerator in 
the displacer (phosphor, tin bronze or stainless steel wire screen), a 
displacer cylinder, a spring and a heat exchanger. The displacer 
with the regenerator was actuated by the pressure difference be¬ 
tween the expansion space and the warm space. In order to reduce 
the mass of displacer, the non-metal material had been adopted to 
process the shell of regenerator. Meanwhile, the axial heat conduc¬ 
tion had almost been avoided due to very low coefficient of heat 
conductivity. Moreover, in order to reduce the axial conductivity 


Compressor Drive Signal 



Linear Compressor Module Expander Module 


Fig. 2. Schematic diagram of Stirling cryocooler. 
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loss of the displacer cylinder, the high strength metal thin wall 
tube with very small section area had also been used. In terms of 
the regenerator, material property, mesh, porosity, weave style 
were the main factors should be considered. Pressure drop and 
heat transfer characteristic of regenerator related closely with 
above parameters could affect the performance of cryocooler obvi¬ 
ously. Determination of these parameters highly depended on the 
accumulated data and empirical correlation expressions [8]. 

One of key components to the overall efficiency of a cryocooler 
was how well heat was transferred from the heat load to the he¬ 
lium in the cold end of cooler. The more efficiently this heat trans¬ 
fer occurred, the lower the temperature difference between the 
heat load and helium in the cold end would be. This resulted in 
the helium not having to be as cold as more efficient cooler. There 
were essentially three resistances to this flow of heat: the first one 
was the interface resistance between the heat load and the heat 
acceptor portion of the cold end, the second one was the conduc¬ 
tion resistance across the heat acceptor, and the last one was the 
convection resistance from the heat acceptor to the helium gas in¬ 
side the cold end. In cryocooler, the first resistance was mainly 
composed of the contact thermal resistances between the heat load 
and cold end, which could be reduced by coating the heat transfer 
enhance grease or cryogenics glue. The red copper cold end with 
very high coefficient of thermal conductivity was welded directly 
with the displacer cylinder. It was the last resistance that is the 
dominant term in the overall resistance, this resistance was de¬ 
fined as 1/(M), where h = convective heat transfer coefficient, 
and A = heat transfer area. A special small-hole structure was 
adopted to increase the velocity of helium scouring the inner face 
of cold end. So the convective heat transfer was enhanced. 

The reliable, lightweight, compact, high efficiency electronic 
controller installed in a separate box contains a moving mass posi¬ 
tion signal processor for the compressor and displacer respectively, 
power supplies, telemetry interface, vibration cancellation cir¬ 
cuitry, and logic circuit. There are two control models could be 
chose. The first one was the close loop control which could be used 
in practical engineering. After inputted the digital single from A/D 
converters of the two LVDT sensors, the servo control system, in 
which a single-chip microcontroller (8031) had been used as the 
main controller, determine the operating parameters of the moving 
masses including direct current bias, frequency, phase between 
compressor and displacer. These parameters then were processed 
by program stored in the chip, which followed the structured pro¬ 
gramming method and includes start-up block, main control block, 
pressure force process block and self-test block telemetry block. 
After iteration calculation, the feed back signal was generated 
and subsequently input parameters could be regulated. The other 
was open loop control model, which was more suitable to the 
experimental study. Many above feedback processes and self-reg¬ 
ulating functions had been cancelled. Various parameters could 
be adjusted manually [9]. 


3. Experimental facility and procedure 

3.1. Experimental facility 

The tested machine is a one-stage Oxford split-Stirling cryoco¬ 
oler constructed at Cryogenic Laboratory in Shanghai Institute of 
Technical Physics (SITP), whose geometrical dimensions and 
operation conditions are list in Table 1. Fig. 3 shows a schematic 
diagram of the experimental apparatus of the FPFD Stirling cryoco¬ 
oler. Two piezoelectric dynamic pressure sensors were used to 
monitor the pressure oscillations at the outlet of the compressor 
and buffers of displacer. Meanwhile, two linear variable-differen¬ 
tial position transformers (LVDT) were mounted at the end of the 
compressor and displacer respectively for the displacement mea¬ 
surement of the pistons, enabling the volume change to be deter¬ 
mined subsequently according to the driven mechanism. A 
platinum resistance thermometer being calibrated was attached 
to the cold head to measure the temperature of the cold end. Film 
heater attached on the cold end has counteracted the refrigeration 
capacity of cryocooler that is equal to the power of heater. The 
apparatus of the Stirling cryocooler, except for the components of 
the room temperature region, was connected to the vacuum flange 
during the experiment; a vacuum chamber was connected to a 
high vacuum pump under a pressure of 10~ 5 Pa. The high vacuum 
pump system consists of a rotary roughing pump, a turbo-molecu¬ 
lar pump and vacuum gauges. 

The following operation parameters can be adjusted, including 
voltage amplitude, voltage frequency, charged pressure, tempera¬ 
ture of cold tip and refrigeration capacity. Moreover, the geometri¬ 
cal parameters including the mesh of regenerator, length of 
connecting tube and the stiffness of spring can also be adjusted 
according to the experiment demands. In this experiment study, 
in order to investigate the characteristic of phase shift further, be¬ 
sides the above common measuring and regulating measures, the 
following special methods have been adopted. 

(1) An extra pneumatic damping ring has been designed and 
installed on the expander assembly. Different damping coef¬ 
ficient can be obtained by adjusting the width of clearance 
cut on the ring. 

(2) In order to get the different natural frequency of expander 
assembly, Additional mass has been screwed on to the 
expander assembly. 

The following tests were undertaken for the experimental per¬ 
formance analysis and comparative study of cryocooler. 

(1) The amplitudes of the displacer’s stroke and phase shift 
under different temperatures. 

(2) The cooling capacity characteristics under the natural fre¬ 
quency of the displacer. 


Table 1 

Dimension and operating conditions of FPFD Stirling cryocooler. 


Working fluid 

Charge pressure (p) 

Working conditions 

Helium 

0.8-1.5 MPa 

Ambient temperature 298 K, relative humidity 60%, cold tip temperature 80 K~200 K 


Parameters of compressor 


Parameters of displacer 


Dimensions of compressor 

0 120 mm x 170 mm 

Dimensions of displacer 

0 160 mm x 148 mm 

Stiffness of piston spring (k p ) 

1.8 N/mm 

Stiffness of displacer spring (I< d ) 

3 N/mm 

Mass of piston (m p ) 

120 g 

Mass of displacer (m d ) 

22 g 

Input voltage 

28 V 

External Diameter of connecting tube 

2.4 mm 

Input excitation frequency (f) 

Specific thrust 

35-44 Hz 

14N/A 

Length of connecting tube 

220 mm 
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Fig. 3. Experimental apparatus of the FPFD Stirling cryocooler. 


(3) The displacer’s amplitudes under various charging pres¬ 
sures, the amplitudes of pressure and of operating 
frequency. 

(4) The phase shifts between the displacements of the piston 
and the displacer under different operating frequencies. 

3.2. Experimental procedure 

In the process of experiment, change of charging pressure and 
additional mass of displacer involve assembly and disassembly of 
cryocooler. In reality, mechanical damping is very sensitive to 
the assembly technology. So in order to guarantee recurrence of 
experiment, the damping characteristic of piston and displacer 
must be measured and subsequently adjusted under the unload 
conditions. Indispensable steps should be implemented in follow¬ 
ing sequence, such as assemble the Stirling cryocooler, damping 
test, leak detection, evacuating cryocooler, gas charging, evacuat¬ 
ing vacuum chamber, connecting the cryocooler with the measure 
instrument and control box, regulating the operating parameters, 
collecting the performance parameters, finishing the experiment, 
closing the electrical power, rewarming cryocooler, change 
components. 


Moreover, this simple description of the pressures is only valid if 
all constants I< p , I< d , C pc , C dc , C pe and C de are independent of x d 
and x p , respectively. Actually, from our previous experimental data 
and published information [11], the stiffness of gas spring (K p ) in 
compression space is so big that has same magnitude with stiffness 
provided by mechanical spring. In addition, this stiffness is very 
sensitive to the stroke of piston. In order to reflect this character¬ 
istic in Eq. (8), the K p should be divided into two parts in reference 
to Eq. (3), the first one (K p = P mean A 2 p /V 0 B) has no relationship with 


the piston stroke; another one (K p = P m A 2 


_ cc _ 

71p(X m —Xp)+V 


1 

V 0 B 


) is the 


function of piston stroke. In contrast, the gas spring effect is far 
weaker in displacer than that in piston, so K p in Eq. (9) equal 
approximately to I< p . So subtracting P c from P e , the pressure differ¬ 
ence is 


Pe Pc — (^pe Cp C )Xp + (C de C dc )X d Xp(X p )X p 

— C pp x p + C dd x d — /C p (X p )x p (10) 

The ideal refrigeration effect and ideal power input, the expressions 
for which are given by Jonge [3] can be rewritten as 

Qeo = n <f> p d dv d = -7rf/C p X p X d A d sin iA pd (11) 


4. Theory of phase shifts 

4.1. Displacement phase shift between displacer and piston 

Motions of displacer and piston can approximately be consid¬ 
ered as the harmonic wave, and are given, respectively, by 


X d =X d COS (COt + l/'pd) 

(6) 

Xp = X p cos (cot) 

(7) 

According to the theory of free piston Stirling refrigerator developed 
by Jonge [10], the pressure fluctuation in the variable compression 

and expansion volumes are expressed as 

Pc = Pm + X p x p + XdX d + C pc x p + C dc x d 

(8) 

Pe = Pm + XpXp + X d x d + C pe Xp + C de X d 

(9) 


This assumption is based on the ignorance of the effect of leak¬ 
age, temperature variations on the pressure. Also the flow is re¬ 
garded as the laminar flow so as to be liberalized conveniently. 


Pd = n j p p dVp = nfK d X p X d A p sin ^ pd (12) 

The ideal coefficient of performance (COP idea i) may now be written 
as 

COPideai = Qeo/Pci = K p A d /K d A p (13) 

It should be noted that, though the refrigeration effect is directly 
proportional to sim /r pd , the ideal COP idea 1 is independent of \j/ pd . 

The area of displacer rod is comparatively smaller than that of 
displacer, so if the area of displacer rod (A r ) is ignored, the Eq. (1) 
can be transform into 

M d X d = A e (p h — P e ) — ^md^d — C md X d (14) 

in which P h can be expressed as P h = P c - A p hc . As far as the con¬ 
struction is concerned, the pressure difference (zlp hc ) between hot 
and compression chamber can be attributed to two factors, the first 
one is the friction loss in the connecting tube, the second one is the 
compressibility and flow loss effect of gas in dead volume if it has. 
According to the published investigation [11] the pressure ampli- 























































































































X. Chen et al. / Cryogenics 49 (2009) 120-132 


125 


tude declines about 5% across the connecting tube, which can al¬ 
most be ignored in this calculation. So the Eq. (14) can be transform 
into 


M d X d — A e (p c P e ) kmdXd Qnd^d (15) 

Substitute Eq. (10) into Eq. (15), the expression is 

M d X d H - ^d-^e^d H - ^md^d — -^e^ppXp T - A e Kp (X p )X p (15) 


in which C d = C md /A e + C dd . Avail of the harmonic wave assumption 
(Eqs. (6) and (7)), the displacement phase shift and amplitude ratio 
can be deduced from the Eq. (16) which is as follow: 


i// pd = arctan 


f (co 2 d - co 2 )M d \ 
v coC d A e ) 


+ arctan 





Mjj(0 - CO 2 ) 2 + (ft)C d /ie) 2 
(COCppAe ) 2 + (M e ) 2 


1/2 



From previous viewpoint [10], only under the conditions of co < co d 
and \j/ pd > 0 can the machine run as a refrigerator. But according to 
experimental phenomenon, the output refrigeration capacity is 
actually not decline to zero, which contradicted with above theoret¬ 
ical prediction. The second term of right side of Eq. (17) explains 
this controversy, which leads to a positive angle and has a positive 
relationship with the stiffness of gas. 


\l/ v „ = arctan 


pp 


t^5" mp X p F e sin y 
_(^mp co 2 M p )X p F e cos y_ 


\F\ 


Xp k mp 



MpCP 2 

^mp 


F e cos y 

^mpX p 


+ 


coC mp F e sin y 


k 


mp 


^mp^p 



( 21 ) 


= fPcdVc = n(coC mp Xp - F e siny) 

PmVm^) 2 (/<mp-C0 2 M p )Xp-F e COSy 

From above equations, the relationship ofi/^ pp and D is D = 7itgil/ pp 
the relationship of ^ pp and F is F = seci^ pp . These expressions reveal 
that the dimensionless amplitude and dissipation is the monodr- 
omy function of phase shift. With the increase of the i/r pp , the 
non-dimension damping (D) will increase correspondingly, mean¬ 
while the non-dimension amplitude (F) will decrease. 


4.3. Relationships between the displacer displacement and system 
pressure 

Eq. (14) can also be used to investigate the phase shift between 
the system pressure and displacer displacement. 

H - k md x d + C md x d = A e (p h — p e ) (23) 


4.2. Phase shift between the system pressure wave and displacement 
of piston 

With respect to the structure of Stirling cryocooler, the pressure 
fluctuation attributed to four predominant causes [10], there are 
piston position determining the total volume, displacer position 
determining the amount of gas in the expansion and compression 
spaces, piston and displacer movement giving rise to gas flow, 
causing pressure differences in the space. Similarly, the pressure 
wave shift can also attribute to four principle reasons, such as 
property change with the temperature, energy hysteresis and 
damping effect as well as energy store and conversion effect. 

In gas spring model [2], gas spring force and displacement are 
related to pressure and volume .the differential equation that is 
transformed from Eq. (2) describing the performance of the mod¬ 
elled gas spring more clearly is 

F p = k m pXp + C mp x p + M p x p — fe (19) 

In which f e = F e cos(cot + y).Substituting the sinusoidal displace¬ 
ment of piston into above equation, the gas spring phase, magni¬ 
tude and dissipation are thus found to be 


From Eq. (10), pressure difference (P h - P e ) is actually the harmonic 
function of displacement of piston and displacer as well as their first 
order differential. Appling the Fourier expansion, this pressure dif¬ 
ference can approximately expressed as 

Ph - Pe = p he sin cot (24) 


Substituted the Eq. (24), the phase shift is 


i^ dd = arctan 


k 


md 


CO 2 Mr 


coC 


= arctan 


M d (( 0 2 - co 2 ) 


md 


ojC 


md 


(25) 


5. Experimental results and discussion 

Fig. 4 shows the real displacement wave of piston and displacer. 
According the characteristic of moving coil linear motor, the excit¬ 
ing force (Bil) has positive relationship with the input voltage. So 
by adjusting the output voltage of regulated power source, the 
amplitude of compressing piston can be controlled. With respect 
to the frequency of displacement, because the Stirling cryocooler 
can actually be regarded as the single input and two output vibra¬ 
tion system, the frequency of response (x d ,x p ) equals to that of 
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Fig. 4. Real time displacement wave of the piston and displacer, when (A) Operating frequency is 41.06 Hz, width of damping clearance is 300um, 10 bar (B) Operating 
frequency is 39.08 Hz, width of damping clearance is300 pm, 14 bar. 
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A B 



Fig. 5. (A) Real wave of displacement of piston and gas pressure, when the operating frequency is 37 Hz, width of clearance is 450 pm. (B) Force vector diagram of piston. 


exciting force (Bil). However due to the strong nonlinear pneu¬ 
matic force driving the movement of displacer, the distortion 
caused by high frequency harmonic wave has happened clearly. 
The phase shift of these two displacement influenced by the geo¬ 
metrical and operating parameters can be measured easily form 
above figures. 

Fig. 5A shows the displacement wave of piston lagging behind 
the gas pressure wave. Similar to the displacement wave of piston 
and displacer, the harmonic deformation of pressure exists. De¬ 
duced from Eq. (2), the vector diagram of force applied on the pis¬ 
ton (shown in Fig. 5B) reveals the relationship of magnitude and 
phase angle of them. It is clear that the inertia force parallel with 
the spring force that is vertical to the damping force. The phase 
shift {il/pp) in Fig. 5B can, respectively influence the spatial position 
and magnitude of force vectors. From this oscillogram, when the 
operating frequency is 37 Hz and width of clearance is 450 um, 
the phase shift is 46.6°. 

5.2. Effect of cold tip temperature on the phase shift 

In the cooling process, with the decrease of the cold tip temper¬ 
ature, average temperature of helium in the regenerator decline 
gradually, which result in the declination of pneumatic driving 
force determined by the pressure drop across the regenerator. As 
a consequence, the stroke of displacer decreases consequently. 
For instance, under the ambient temperature (293 K) dynamic vis¬ 
cosity of helium is 19.52 x 10~ 6 N S/m 2 and the density is 2.377 kg/ 
m 3 . In contrast, when the temperature is 80 K, dynamic viscosity of 
helium become 15.2 x 10 -6 N S/m 2 , which decrease about 22% and 
density become 8.936 kg/m 3 , which is 3.759 times of that of ambi¬ 
ent temperature. When the temperature of cold tip decreasing 
from 230 to 63 I< and keeping constant piston stroke (3.8 mm), as 
a result of variation of property displacer stroke decrease corre¬ 
spondingly from 1.45 to 0.95 mm. So in order to shorten the period 
of cooling process in practical operation, the piston stroke should 
be increased by request of enough displacer strokes. 

In terms of displacement phase shift between the piston and 
displacer, due to spring stiffness sensitive to piston stroke, natural 
frequency coupling with damping coefficient (C pp ) [11] increasing 
simultaneously with the decrease of cold tip temperature, in other 
words, numerator and denominator of the second term in Eq (17) 
increasing at the same time, the phase shift just increase in a little 
degree consequently. For example, for working condition (A), when 
temperature decrease from 270 to 60 K, the phase shift vary from 


87.5° to 89.5°, just increase 2.28%. From Eqs. (11) and (12), if ignor¬ 
ing the effect of varying phase shift, in cooling process the refriger¬ 
ating capacity and efficiency has positive relationship with the 
magnitude of stroke approximately. So for Striling cryocooler, if 
the low temperature were needed, the enough strokes of piston 
and displacer should take measure to be provided and the collision 
between piston and cylinder should be avoided simultaneously. 

The resonance characteristics [12] of the compressor also 
change somewhat between ambient and cryogenic temperatures, 
which meanly attributes to the stiffness of the gas spring drops 
with decreasing gas pressure as gas density increases in the cold 
finger. Note that the damping contrastingly increases somewhat 
at cryogenic temperatures. Meanwhile the damping also increases 
somewhat with increasing stroke. So as shown in Fig. 6, the phase 
shift increase with the decrease of temperature following the prin¬ 
ciple given by Eq. (20). Furthermore, the non-dimension damping 
coefficient given by Eq. (22) increases about 10%. So it can be con¬ 
cluded that the working performance actually deteriorate under 
the low temperature working condition. 

With reference to the displacer phase, C md will increase in con¬ 
sequence of lowing temperature. Meanwhile the gas spring stiff¬ 
ness of displacer keeps almost constant, so the phase will 
decrease as shown Fig. 7. The increasing ratio of damping can be 
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Fig. 6. Displacement phase shift between displacer and piston & phase shift 
between system pressure and the displacement of piston with the variation of cold 
tip temperature. 
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Fig. 7. Phase shift between displacer displacement and the system pressure with 
the variation of cold tip temperature. 

deduced from the variation of phase shift. For instance, for condi¬ 
tion (A), when temperature varying from 270 to 60 K, the damping 
coefficient (C md ) increase 1.6885 times. 

From another viewpoint, it can be found that as for three phase 
shifts, that of the displacer with wider clearance (450 pm) is bigger 
than that with 300 pm clearance. It reveals that the less the system 
damping has, the bigger phase shift will be. 

In general, cold tip temperature causing the variation of helium 
property increases damping of both piston and displacer. But it 
seems that the damping of piston is more sensitive to the stroke 
than property variation. And movement characteristic of piston is 
controlled in larger degree by the external factors including elec¬ 
tric voltage and frequency. In contrast, the motion of displacer de¬ 
pends clearly on damping coefficient instead of stroke. 

5.2. Effect of operating frequency on the phase shift 

Frequency adjustment can be realized conveniently by the var¬ 
iable frequency power source. In this series experiment, only fre¬ 
quency of input electric power is adjusted. For the displacer with 
300 pm clearance, while frequency varying from 43 to 37 Hz, the 
strokes of piston increases from 2.36 to 3.983 mm correspondingly. 
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Fig. 8. Displacement phase shift between displacer and piston and phase shift 
between the system pressure and displacement of piston with the variation of 
working frequency. 


So referring to the displacement phase shift, although the damping 
(C d ) decrease with the declination of stroke of displacer, working 
frequency is actually the main influencing factor which can be de¬ 
duced from Eq. (17). As shown in Fig. 8, for the displacer with 
450 pm clearance and 46 Hz natural frequency stroke phase de¬ 
crease from 115° to 92.6° while working frequency varying from 
37 to 43 Hz. When the working frequency is 41 Hz, it also can be 
drawn from Fig. 8 that the more natural frequency and wider clear¬ 
ance has, the bigger phase shift will be. This tendency obviously 
coincides with that predicted by Eq. (17). 

Conversely, the pressure phase shift increase slightly with in¬ 
crease of working frequency. For instance, the phase shift increases 
about 5.15% when working frequency increase from 37 to 43 Hz for 
displacer with 450um clearance and 46 Hz natural frequency. For 
compressor, increasing frequency brings on decreasing stroke (X a ), 
subsequently the damp (C mp ) and stiffness of piston (J( mp ) will de¬ 
cline. From Eq. (20), the comprehensive effect of above factors leads 
to increase of phase shift. Meanwhile, the almost fixed the dimen- 
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Fig. 9. Phase shift between the displacer displacement and system pressure with 
the variation of working frequency. 
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Fig. 10. Displacement phase shift between displacer and piston and phase shift 
between the system pressure and displacement of piston with the variation of 
charge pressure, while f n = 46 Hz, W c \ = 300 pm. 
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sionless damp shown in Eq. (22) accounts for the performance of 
compressor just decline a little with variation of working frequency. 

From classical vibration theory, when working frequency is less 
than natural frequency, the displacer phase shift will rise with in¬ 
crease of working frequency. Higher natural frequency or lower 
damping coefficient will lead to bigger phase shift. It also can be 
seen that the resonant frequency is actually less than the natural 
frequency and the difference between them depends on damping 
coefficient. From Fig. 9, it is clearly that less damping and large nat¬ 
ural frequency leads to less phase shift, which conforms to the pre¬ 
diction of Eq. (25). However, it exists a paradox seemingly between 
Fig. 9 and Eq. (25) that when the working frequency is 43 Hz, the 
phase shifts of displacer whose natural frequency is 43 Hz do not 
equal to zero (Fig. 9), which can be deduced from Eq. (25). This con¬ 
tradiction actually attribute to the traditional calculation method 



Charge Pressure(bar) 


Fig. 11. Phase shift between the displacer displacement and system pressure with 
the variation of charge pressure, while f n = 46 Hz, W d = 300 pm. 




used in primary design, that is for displacer the mechanical spring 
stiffness is just involved in the natural frequency calculation. But 
in fact, the gas in the hot and expansion chamber can provide stiff¬ 
ness to some extent, which is far less than that of mechanical spring 
and ignored at most cases. This gas stiffness is also the function of 
working parameters. 

5.3. Effect of the charge pressure on the phase shift 

For displacer with f n = 46 Hz, W d = 300 pm, when working fre¬ 
quency is 41 Hz and charge pressure varies from 9 to 14 bar, the 
stroke of piston decrease correspondingly from 3.93 to 3.6 mm 
and the stroke of displacer also decrease about 15%. This tendency 
can also applied to other working frequencies [12,13]. From Eq. (3), 
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Fig. 13. Displacement phase shift between displacer and piston & phase shift 
between the system pressure and displacement of piston with the variation of 
displacer’s natural frequency, while P c = 15 bar, W d = 300 pm. 
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Fig. 12. When/ n = 46 Hz, W d = 300 pm, with the variation of charge pressure (A) the lowest reached temperature, (B) cooling capacity when the temperature is 80 K, (C) input 
power when the temperature is 80 K, working frequency is 43 Hz, (D) COP when the temperature is 80 K, working frequency is 43 Hz. 
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the stroke and charge pressure has opposite effect on the gas stiff¬ 
ness and damping coefficient. However as shown in Fig. 10, these 
two phase shift have the positive relationship with the charge 
pressure. Pressure phase shift shown in Eq. (20) have negative rela¬ 
tion with the piston stroke. As for the gas stiffness of piston, the 
increasing mean pressure and dead volume have conversely effect 
on the value of stiffness (seen in Eq. (3)). So the gas stiffness of 
compressor keeps constant as a whole. On the contrary, natural 
frequency of displacer has not same sensitivity to stroke as that 
of piston. From Figs. 10 and 11, it can be deduced that increasing 
phase shift imply that the effect of decreasing stroke of displacer 
on the damping coefficient is bigger than that of increasing mean 
pressure. In the end, the displacer damping will decline. For in¬ 
stance, from Fig. 11 when working frequency is 41 Hz, with the 
variation of charge pressure from 9 bar to 15 bar, the displacer 
phase shift increase from 38.3° to 54.4°. The corresponding damp¬ 
ing decrease ratio is 56.54%. 

The effect of charge pressure on performance is shown in 
Fig. 12A-D, respectively. By means of adjusting working parame¬ 
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Fig. 14. Phase shift between the displacer displacement and system pressure with 
the variation of natural frequency of displacer, while P c = 15 bar, W c \ = 300 pm. 


ters repeatedly, the lowest temperature obtained under 15 bar is 
51 K. When increasing charge pressure, the amplitude of pressure 
wave will also increase which leads to the increase of cooling 
capacity and input power simultaneously. In reference to Figs. 10 
and 11, the increasing phase shift also contributes to the rise of 
cooling capacity and input power. But from Eqs. (11) and (12) 
the COP is just the uniform function of charge pressure. So charg¬ 
ing high-pressure gas is an effective way to get high efficiency if 
the ability of pressure resistance is permitted, which depend on 
the thickness of shell and material property. But it should be point 
out that increasing rate of COP slow down when the charge pres¬ 
sure is high (as shown in Fig. 12D). The evaluation of advantage 
and disadvantage of increasing of charger pressure is naturally 
very necessary for an optimum design. 

5.4. Effect of the natural frequency on the phase shift 

In this study, in order to change the natural frequency of dis¬ 
placer, the additional mass was attached to the displacer. 1.95 g 
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Fig. 16. Displacement phase shift between displacer and piston and phase shift 
between the system pressure and displacement of piston with the variation of 
clearance width, while P c = 15 bar, f n = 46 Hz. 




Fig. 15. When P c = 15 bar, W d = 300 pm, with the variation of natural frequency (A) the lowest reached temperature, (B) cooling capacity when the temperature is 80 K, (C) 
input power when the temperature is 80 K, working frequency is 43 Hz, (D) COP when the temperature is 80 K, working frequency is 43 Hz. 
















130 


X. Chen et al./Cryogenics 49 (2009) 120-132 


and 7.76 g additional weight parallels to 43 Hz and 40 Hz natural 
frequency, respectively. 

As shown in Fig. 13, in terms of the displacement phase shift, for 
different three natural frequency of displacer the higher the working 
frequency is, the lower phase shift will be. If co = co d (displacer is 
turned to the running frequency), then cp pd ^ 0, which do not coin¬ 
cide with the prediction given by Jonge [10]. Actually, the theory 
developed by Jonge is the special case of Eq. (17) when the gas spring 
is ignored. In the real test, compressor running full stroke is very 
dangerous, at this moment the autocontrol system and anti-colli¬ 
sion function has been cancelled and substituted by manual opera¬ 
tion. So in order to avoid the damage of collision, the full stroke is 
practically rare. In this case, when the running frequency equal to 
the natural frequency, positive displacement phase shift still exist 
and have inverse proportion with the working frequency. For in¬ 
stance, when natural frequency of displacer is 43 Hz, displacement 
phase shift is 35.5° under the condition of 43 Hz working frequency. 

Fig. 13 indicates that pressure phase is not sensitive to the 
working frequency. And it can also found that when natural fre¬ 
quency of displacer is 43 Hz, the maximum pressure phase shift 
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Fig. 17. Phase shift between the displacer displacement and system pressure with 
the variation of clearance width, while P c = 15 bar, f n = 46 Hz. 


(44°) happens under the condition of 41 Hz working frequency. 
From Eq. (20), when working frequency of compressor equals to 
that of displacer, the pressure phase reach its peak. According to 
the vibration theory, the resonant frequency is actually less than 
natural frequency, which implies natural frequency of piston less 
than 41 Hz. As a result, it can be drawn from Fig. 13 that natural 
frequency of compressor should be less than that of displacer, 
which is the basic principle followed in real design (Fig. 14). 

As far as expansion phase is concerned, when natural frequency 
of displacer is 40, 43 and 46 Hz respectively, 37 Hz working fre¬ 
quency can lead to comparative bigger phase shift to different nat¬ 
ural frequencies than 39, 41 and 43 Hz. In this case, the frequency 
ratio (cDolcDn) is 0.925, 0.86, 0.8 correspondingly and 0.86 can bring 
up maximum phase shift that approach to the optimum range of 
frequency ratio (0.75-0.85) provided by Park [14]. Furthermore, 
the reason why expansion phase shift do not become zero when 
working frequency equal to resonant frequency is gas stiffness of 
expansion, especial hot space do not involve into the calculation 
of natural frequency. Similar to displacement phase shift, this igno¬ 
rance can cause additional phase shift. 

Fig. 15 shows that natural frequency has obvious effect on the 
performance of cryocooler. For the displacer with natural fre- 
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Fig. 18. When P c = 15 bar,/ n = 46 Hz, with the variation of clearance width (A) the lowest reached temperature, (B) cooling capacity when the temperature is 80 K, (C) input 
power when the temperature is 80 K, working frequency is 43 Hz, (D) COP when the temperature is 80 K, working frequency is 43 Hz. 
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Fig. 20. (A) Cooling temperature vs. cooling power and (B) cooling temperature vs. COP. 


quency 46 Hz, when working frequency is 43 Hz, the lowest 
reached temperature is 51 K and cooling power and input power 
at 80 K is 0.64 and 27.25 W, respectively. At this point, the dis¬ 
placement phase shift is 71.3°. In refer to Eqs. ((11)—(13)), the in¬ 
crease of amplitude of piston (X p ) and displacer (X d ) lead to 
decrease of displacement phase shift ( 0 ), similarly the increase of 
pressure wave causes the decrease of amplitude. So if the maxi¬ 
mum cooling power is wanted, the working frequency should be 
close to the natural frequency of displacer. 

5.5. Effect of the damping on the phase shift 

Fig. 16 indicates that when working frequency is 41 Hz, dis¬ 
placement phase shift varies from 85.2° to 101.8° with the increase 
of clearance width from 200 to 450 pm. As predicted by Eq. (17), 
increasing damping of displacer cause decreasing displacement 
phase shift. As for pressure phase shift, when damping of displacer 
increases as a result of decreasing clearance width, from Eq. (18), 
the amplitude ration will decrease; the stroke of piston and dis¬ 
placer will all increase consequently. In the end, according to Eq. 
(20), the pressure phase shift will increase as shown in Fig. 16. 
Experimental data finally verified this conclusion. For instance, 
when working frequency is 41 Hz and width of clearance vary from 
200 to 450 pm, the piston stroke increase from 3.76 to 3.84 mm, 
correspondently the stroke of displacer increase 1.923 times. As 
for expansion phase shift, with the increase of clearance, the damp¬ 
ing provided decrease; in consequence the phase shift increase. 
Discrepancy between experiment data and theory attributes to 
the rigging error (see Figs. 17 and 18). 

For 43 Hz working frequency and 46 Hz natural frequency, the 
lowest temperature obtained is 51 K. Cooling power and input 
power for 80 I< cold tip temperature is 0.635, 27.2 W respectively. 
From Eqs. (11) and (12), the cooling power and input power have 
actually close relationship with displacement of piston, displacer 
and phase shift between them, all of which can be influenced by 
the damping coefficient of displacer. Maximum COP 2.3 indicates 
that 300um clearance is optimum width in comparison to others. 

5.6. Optimum operating condition of Stirling cryocooler 

From above experimental analysis, the optimum working con¬ 
ditions for 80 K cold tip temperature is charge pressure 15 bar, nat¬ 
ural frequency of displacer 46 Hz, width of clearance 300 pm and 
working frequency 43 Hz. At this condition 0.632 W cooling power 


is obtained at the cost of 27 W input power and COP is 2.31%, low¬ 
est obtained temperature is 51 K. Fig. 19 shows cooling curve verse 
time under optimum operating conditions. It can be seen that from 
ambient temperature to 77 K, 22 min is needed; 50 min is spent to 
get the lowest temperature 51 K. The cooling rate of initial stage is 
faster than that of later stage. Fig. 20A shows the effect of cooling 
temperature on cooling power under the 26 W input power, the 
approximate linear relationship exist between them and slope ra¬ 
tio of this line is 0.02 W/I< within the range from 71 to 100 K. Sim¬ 
ilar to cooling power, COP also is almost linear function of cooling 
temperature at same temperature range. The decline rate of COP 
against cooling temperature is 0.0794%/K. 

6. Conclusions 

Phase shift characteristics including displacement, pressure and 
expansion phase difference have been studied theoretically and 
experimentally, respectively. Here, the main conclusions have been 
drawn in the following: 

(1) According to demands of maximum cooling power and COP, 
the optimum working conditions of this Stirling cryocooler 
for 80 K cold tip temperature are as follows: charge pressure 
15 bar, natural frequency of displacer 46 Hz, width of clear¬ 
ance 300 pm and working frequency 43 Hz. Compared with 
working frequency, stroke phase and some other parame¬ 
ters, actually phase shift is the most complicated one involv¬ 
ing relative movement of two pistons. So in the process of 
optimum seeking, an indispensable step is to locate the 
working point in the optimum range of phase shift that 
can actually act as the index of optimization. As for displace¬ 
ment phase shift, neighborhood interval of 90° is the ideal 
working domain. Meanwhile, the expansion phase shift 
should be adjusted approaching to 0° as near as possible. 
Pressure phase shift, which reflects the damping characteris¬ 
tic of compressor, should also be as small as possible. 

(2) Eqs. (17), (20) and (25) reveal the relationship between three 
phase shifts and working parameters respectively, predica¬ 
tion from which agree with the experimental result qualita¬ 
tively. Due to the difficulties of measuring pneumatic 
damping coefficient precisely and coupling characteristic of 
different parameters highly, quantitative analysis must 
resort to abstruse mathematic method, which is very hard 
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to be used in practical adjustment. Conversely, these equa¬ 
tions describe the phase shifts function clearly and are able 
to guide the optimum seeking directly and effectively. Fur¬ 
thermore, the discrepancy between previous theory and cur¬ 
rent experimental phenomenon has been avoided in Eq. (17) 
by involving the gas spring stiffness into the mathematic 
model. The dimensionless damping has been expressed as 
the function of pressure phase shift, which can be served 
as the index of performance of compressor. 

(3) By experimental study, many running characteristics have 
been inferred from the variation of phase shift. The variation 
of helium property caused by decreasing cold tip tempera¬ 
ture increases the damping of both piston and displacer. 
But it seems that the damping of piston is more sensitive 
to stroke than property variation. And movement character¬ 
istic of piston is controlled in larger degree by the external 
factors including electric voltage and frequency. In contrast, 
the motion of displacer depends clearly on pneumatic damp¬ 
ing coefficient instead of stroke; the increasing working fre¬ 
quency results in stroke decrease of piston and displacer. In 
terms of piston, the effect of increasing frequency counter¬ 
acts decreasing damping caused by the decreasing stroke. 
In contrast the movement of displacer is influenced by fre¬ 
quency to larger extent. As the charge pressure is increased, 
resonant frequency of the compressor is accordingly 
increased on the condition of constant piston stroke, but 
the resonant frequency of the displacer, which is decided 
by the mechanical spring, is almost constant. And natural 
frequency of displacer has little effect on the movement of 
piston. If damping coefficient of displacer keeps constant, 
working frequency should match the natural frequency; 
damping of displacer has negligible impact on the piston 
movement. There exists the matching problem between 
the natural frequency and damping coefficient. Higher natu¬ 
ral frequency of displacer should be adopted if the damping 
coefficient is great. 
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